Neel and disordered phases of coupled Heisenberg chains with S — \ to S — 4 
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We use the two-step density-matrix renormalization group method to study the effects of frustra- 
tion in Heisenberg models for 5*=! to .5 = 4 in a two-dimensional anisotropic lattice. We find that 
as in S = | studied previously, the system is made of nearly disconnected chains at the maximally 
frustrated point, Jd/J± = 0.5, i.e., the transverse spin-spin correlations decay exponentially. This 
leads to the following consequences: (i) all half-integer spins systems are gapless, behaving like a 
sliding Luttinger liquid as in S = 4; (ii) for integer spins, there is an intermediate disordered phase 
with a spin gap, with the width of the disordered state is roughly proportional to the ID Haldane 
gap. 



I. INTRODUCTION 

An important theorem by Dyson, Lieb, and Simon 
(DLS) states that the Heisenberg Hamiltonian on bipar- 
tite lattices with S > 1 have long-range order in the 
ground stated. We now know from quantum Monte 
Carlo simulations^ that this result extends to S = | sys- 
tems. There is current interest as to how a disordered 
state emerges out of the Neel state. This question is im- 
portant to the physics of frustrated materials where some 
systems exhibit no magnetic order down to the experi- 
mentally accessible temperatures and are thus believed 
to be disordered in the ground state. A disordered phase 
may also be relevant in the theory of high temperature 
superconductors. This issue and the related one of the 
eventual role of Berry phases was hotly debated in the 
late 1980's and early 1990's soon after the discovery of 
the high Tc materials^SiL&. Many possible disordered 
states were proposed but none of them gained consensus. 
For an extensive discussion on this topic, we refer the 
reader to the book by Fradkin— 

In the search of the nature of a ground state of a 
quantum Hamiltonian, there is another useful theorem 
by Lieb, Schultz, and Mattis (LSM) 2 ^, initially formu- 
lated for ID systems and later extended to 2D systems 
by Affleck 8 restricts the possible ground states of half- 
integer spin systems. Either the ground state is degen- 
erate, (presumably due to a broken symmetry) or it is 
unique and gapless without any long-range order. For 
Heisenberg Hamiltonians, the latter possibility appears 
difficult to realize in 2D because dimer-dimer or spin-spin 
correlations have a power law decay in ID. It would thus 
be expected that interchain couplings will lead to long- 
range order in one of the two channels. If the magnetic 
order is frustrated, a dimerization would be expected. 
Frustrated spin models often involve competitions be- 
tween two magnetic orders, the expected dimerized phase 
is believed to lie between these magnetic phases. A possi- 
ble alternative to this scenario allowed by the LSM the- 
orem is the occurence of a disordered gapless state at 
transition between the magnetic phases. 

In recent publication s 24 ! 25 , we have studied the possi- 
ble emergence of a disordered state in a model of coupled 



Heisenberg chains. This model is a spatially anisotropic 
version of the well studied J\ — J2 model. The model 
essentially retains the physics of the J% — J2 model: it 
presents a phase transition between two magnetic phases. 
The first phase is a Neel phase characterized by the or- 
dering wave vector Q = (tt,tt). The order parameter in 
this phase is maximal in absence of the diagonal coupling 
Jd- It decreases as Jd is increased until it vanishes at the 
maximally frustrated point. Beyond this point, another 
Neel state with Q = (tt, 0) becomes the ground state. 
The existence of these two phases was predicted in nu- 
merous studies^ of the isotropic ( J± = 1) version of this 
model. But what has caused the continuing interest into 
this simple model is the question of whether there is an 
intermediate phase between these two magnetic phases. 
Simple physical arguments suggest the existence of such 
a phase because the two states are associated with sub- 
groups of the SU (2) symmetry group of the spin Hamil- 
tonian and of the C4 symmetry group of the square lat- 
tice that neither include each other. From the Landau 
theory we know that a continuous transition from these 
two phases is forbidden. Initial suggestions for the possi- 
ble intermediate phase was the resonating valence bond 
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chiral spin liquid— or 



the spin-Peierls (SP)5"k The SP phase has lately emerged 
as the front runner. It causes the same type of difficulty 
as the direct transition between two magnetic phases. 
Since in that case, the transition is between a magnetic 
state that breaks the spin rotational symmetry but not 
the lattice translational symmetry to a dimerized state 
that breaks the lattice translational symmetry but not 
the spin rotational symmetry. The transition to SP phase 
has triggered interesting proposal of an extension of the 
conventional Landau-Ginzburg- Wilson (LGW) theory of 
second order phase transitional. 

It seems, however, that this theory does not apply to 
the Ji — J2 model, where we did not find any evidence 
of an intermediate phase for S = 4. Numerical data 
suggests a transition, which seems to be of second or- 
der, between the two magnetic phases at the maximally 
frustrated point for both the anisotropic and isotropic 
models^ 2 -. At the transition point, the competing mag- 
netic orders neutralize each other and the system behaves 
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like a collection of loosely bound chains, even if the bare 
interactions are not small. Classically, the ground state 
is degenerate at this point. This degeneracy is lifted 
by quantum fluctuations. In Refi22i2^, we have shown 
that among all the possible clusters, chains offer the best 
compromise between minimizing the energy and avoid- 
ing frustration at the same time. At the maximally frus- 
trated point, the transverse interactions seem to be irrele- 
vant. I.e., up to the largest lattice size studied, transverse 
spin-spin correlations decay exponentially and the longi- 
tudinal correlations revert to those of decoupled chains. 
This disordered state is a singlet and gapless consistent 
with the LSM theorem&Sii. It is reminiscent of a slid- 
ing Luttinger liquid (SLL) found in models of coupled 
fermions chains^iS*ii. It thus appears that the inter- 
mediate region where a disordered phase has long been 
thought to exist is just a critical region. That is the rea- 
son why it has resisted various approaches for nearly two 
decades. At the maximally frustated point, the correla- 
tion functions are ID-like, thus the rotational spin sym- 
metry of the system is restored. In fact such a transition 
between these two magnetic phases already exists in the 
unfrustred model when the transverse exchange parame- 
ter, J±, is varied from positive values to negative values. 
At the point J± = 0, there is a transition from 2D to ID, 
where properties are identical to that of the maximally 
frustrated point. The only difference between the case 
J± = and the maximally frustrated point is the pres- 
ence of irrelevant transverse terms which do not change 
the long distance behavior of the correlation functions as 
shown in Ref^. From this result the LGW theory ap- 
plies if it is assumed that the system's group of symmetry 
at the critical point contains the groups of symmetry of 
the two magnetic phases. 

It is important to study how this interesting physics 
extends to larger spin systems. First, because many 
frustrated systems contain larger spins. Second, because 
there are some interesting predictions from large S ap- 
proaches about the emergence of a disordered phase from 
a Neel phase as function of S. Affleck 8 argued that since 
the LSM theorem does not apply to integer spin sys- 
tems, there might be a distinction between integer and 
half-integer spin systems in 2D as well. Haldane^ dis- 
cussed that in addition to the now well-established dis- 
tinct behavior between half-integer and integer spins in 
one dimension's, there might be a difference between odd 
and even integer spins in two dimensions due to the ef- 
fects of the Berry phase. Read and Sachdev carried out 
a large N analysis of the possible disordered phase as 
function of the value of the large N equivalent of the 
spin. Their results were consistent with Haldane's pre- 
dictions. Three types of disordered states were predicted. 
For half-integer spins, the non-magnetic phase is a SP 
phase which breaks the translational symmetry along the 
two directions of the square lattice. For odd integer spins, 
the non-magnetic state is made of weakly-coupled chains, 
i.e., the translational symmetry is broken along one direc- 
tion only. Finally for even integer spins, the disordered 



state are valence bond solids, like the Affleck-Kennedy- 
Lieb-Tasaki (AKLT) state 33 , i.e, it does not break any 
translational symmetry. 

In the large S approaches the Heisenberg model is 
mapped onto the non-linear sigma model {NLo) with a 
Berry phase term. This mapping is only approximate^^ 
and there can be some subtle differences with the original 
model 3 ^. Indeed, the realization of the Haldane conjec- 
ture in ID shows their power. But, in absence of ex- 
act results for small S, it is impossible to know whether 
their predictions of a disordered phase extends to small 
S. Another potential problem is that the mapping to 
the cr-model assumes the presence of a smooth configu- 
ration of spins. This is true in the weak-coupling regime 
(Neel ordered phase). But this assumption may break 
down in the strong coupling regime (disordered phase). 
In one dimension, the cr-model coupling constant is given 
by g = 2/(y/s(s+l)y/l - Thus for S = ±, 

the equivalence between the two models breaks down at 
J2 = 0.25, i.e., close to the transition to a dimerized state. 
Such a breakdown seems to occur in the J\ — J2 model 
where, as seen above, the large N predictions conflict for 
S = \ with the TSDMRG in the J x - J 2 model. But since 
spin-half integer systems are critical, it could be objected 
that the behavior seen in our numerical studies are due 
to finite size effects. For large enough lattices there might 
be a relevant interaction which can drive the system to 
a SP phase as predicted from large N. Though this sce- 
nario appears to be unlikely, as discussed in Ref^, it 
cannot be completely rejected. In principle, one would 
expect a different behavior for integer spin systems which 
are known to have a spin gap in IDiS. 

In this paper, we applied the TSDMRG to study 
S = 5, 1, §,...,4 systems in the anisotropic 2D Hamil- 
tonian QJ. Our results in the absence of frustration are 
in agreement with the DLS theorem Refill. We find that 
for all S the ground state is ordered in the absence of frus- 
tration. For all S except for S = h, the order parameter 
is large enough so that the extrapolated value are reli- 
able. This result constitutes a non-trivial test of the TS- 
DMRG, since the TSDMRG starts from decoupled chains 
which arc disordered. When the frustration is turned on, 
the general mechanism found for the destruction of the 
Neel phase is the severing of the frustrated bonds in the 
transverse direction, leading to a disordered state with 
the transverse correlations that decay exponentially at 
the critical point as previously found for S = i. How- 
ever, a different conclusion is to be drawn for half-integer 
S and for integer S for which the LSM theorem does not 
apply. All half-integer systems are similar to S = | . The 
disordered state is confined at the critical point, it has 
a SLL character. But for integer S, because of the Hal- 
dane gap Ah in the chain, there is an intermediate phase 
whose width is roughly oc Ajy. 

This paper is organized as follows. In the next section 
we discuss the model and the method. In section (jlllll . 
we present extensive results for 5 = 1 systems. This 
analysis is similar to the one made for spin i systems 
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in Refi^£. In section IjlVfl the results for systems with 
S = 5 to 4 are presented. In section (0, we present our 
conclusions. 

II. MODEL AND METHOD 
A. model 

We apply the TSDMRG 2 ^ to the spatially 
anisotropic Heisenberg Hamiltonian, 

H = J|| 2^ Si.iSj+i,; + Jj_ ^ Sj^Si^+i 

Z,Z i,l 

+Jd / ^(Sj^Sj-i-i^-!-! + Si+i^Sj^-i-i) (1) 

where J|| is the intra-chain exchange parameter and is 
set to 1; J± and Jd are respectively the transverse and 
diagonal interchain exchanges. This model is the object 
of current interes1*i2i±&i£i2&. It is a starting point of un- 
derstanding the Ji — J2 model which is recovered when 
J11 = J± = J\ and Jd — J%- It retains the basic physics 
of the Ji — J2 model and has the advantage that in the 
limit J_i , Jrf <C 1 , well tested ID results can be used to 
initialize a perturbative RG analysis. 

B. The two-step DMRG 

In the TSDMRG, to study a 2D lattice of size L x (L + 
1) ( we will refer to the 2D systems only by their linear 
dimension L), we start by applying the usual ID DMRG 
(mi states are kept) or exact diagonalization (ED) to a 
single chain I of length L to obtain 7712 low lying eigen- 
states and eigenvalues, 4> ni , £n n n-i — 1,2, ...7712, respec- 
tively. Then, we formally write the tensor product of the 
L + 1 chains, 

*[n] = 4>m^n 2 -4n L+ x- (2) 

is an eigenstate of the Hamiltonian with Jj_ = 0. 
and Jd = 0, [n] = [n\, n-2, til+i)- The $[„] consti- 
tute a many-body basis of the truncated Hubert space of 
the tensor product of L + 1 chains. The corresponding 
eigenvalue is: 



state, the TSDMRG is able to reach the ordered state 
without any addition of a term that explicitly breaks the 
symmetry such as a magnetic field. The TSDMRG was 
tested against the quantum Monte Carlo (QMC) method 
in Ref.( 2 - 5 ) and against ED in Ref.(^). The TSDMRG 
is variational, its performance can systematically be im- 
proved by increasing mi and m.2- Key indicators about 
the performance of the TSDMRG are the truncation er- 
ror p\ during the first step, the width SE of the 777,2 states 
kept and the truncation error p2 during the second step. 
In principle, it is necessary that the ratio of 5E over the 
transverse coupling be large for the TSDMRG to yield 
great accuracy. Typically, one must have SE/J± w 10. 
If this condition is fulfilled and the m,2 states are accu- 
rate enough, .ie., p\ is small, the TSDMRG method can 
reach the QMC accuracy. So far, this has been achieved 
only for small couplings J± < 0.1 and lattice sizes of up 
to L = 16 keeping up to m2 = 96. The amount of calcu- 
lations involved remains modest and so far are done on 
a workstation. The accuracy decreases by increasing Jj_ 
leading to less accurate results in the ordered state. But 
when both J± and Jd are turned on, the performance 
of the TSDMRG becomes more complex as we will see 
below. 

In this work, the calculations for S — 1 were performed 
similarly to those for spin S = h systems in Ref. 22 . In 
most cases, ED was applied during the first step. In some 
case, i.e, some runs of L = 10 and all runs of L = 12, a 
single DMRG iteration was used. For instance for L = 10 
and S = 1, when mi = 81 states are kept, a single 
DMRG iteration is necessary to reach desired size. For 
this calculation, we kept up to m2 = 96 states during the 
second TSDMRG step. For the maximum performance 
of the algorithm, it is necessary that J± be in the order 
of the finite size gap of the single chain A l ■ For S = 1 , 
A L -> A H when L -> oo, where A H = 0.4107iS is the 
Haldane gap. A second condition to fulfill is 5l ^S> J± 
where 5l is the width of the retained eigenvalues. As 
noticed in Refi^, the TSDMRG is more accurate in the 
highly frustrated regime. In Fig.Q we show the trunca- 
tion error, p2, when two states are targeted in the second 
step as function of Jd for S — 1 and J± = 0.4. p2 is 
minimal near Jd — 0.22. At this point the system is 
an assembly of nearly disconnected chains; the DMRG is 
thus expected to perform better. 



C. Illustration in the 5 = 1 case 



E[ n ] — £ni + £n 2 + ••• + £n L + 1 (3) 

The 2D Hamiltonian Q is then projected onto this 
truncated basis to yield an effective one dimensional 
Hamiltonian which is studied using the DMRG. 

The TSDMRG is perturbative; but the expansion is 
made onto the smaller term of the Hamiltonian itself not 
on the Green's function or the ground-state wave func- 
tion. We have shown that starting from a disordered 



We now wish to provide a detailed description of a typ- 
ical TSDMRG calculation. For this illustration, S = 1 
and L = 12, i.c, following the convention set above, the 
size of the 2D lattice is 12 x 13. We start the usual 
ID DMRG iteration keeping mi = 243 states, i.c, the 
initial superblock size is L = 12. At this point the 
DMRG is equivalent to ED. At the next iteration, the 
superblock size is L — 14, its total number of states is 
M s = (3 x 243) 2 . The size of the reduced superblock 
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FIG. 1: Truncation error in the second TSDMRG step as 
function of Jd for J± — 0.4, S = 1 and L — 10, ED was 
performed in the first step and in the second step, 7712 = 64. 



(the superblock minus the two single site blocks) size 
is L — 12. During this iteration, spins sectors with 
Sj. = 0, ±1,±2,±3 are targeted. The lowest states in 
each of these sectors for L = 12 have respectively the 
following energies, E = -16.8696, E ±1 = -16.3854, 
E± 3 = -16.2780. The truncation error is 7.3xl(T 7 . The 
reduced superblock is then diagonalized and the m 2 = 64 
lowest lying states are kept. The energy of the highest 
state among these TO2 states is —14.2687 which is lower 
than E ±i = —12.6524, the lowest states of the = ±4 
sectors. For this reason, these sectors was not targeted. 
The operators Si at each site are stored and updated. 

From the energy levels above, we see that the finite 
size spin gap is = 0.4842, which is not very far from 
its value Ah = 0.4107 in the thermodynamic limit. Our 
choice of J± = 0.4 ensures that the chains will effectively 
be coupled at this size. The matrix O whose columns are 
made of the 772,2 vectors <fi n kept is used to express all the 
operators in the truncated reduced superblock basis, 



diagonal in the basis of the 7712 states kept. The effective 
first order Hamiltonian which approximates the original 
2D Hamiltonian is now given by 

1 1 1 

where, 

S; * Sj + i = ^ Sj^Sj^-i-i (7) 

i 

and 

S; x Sj + i = ^ SuSi+i^+i + Si+i^Su+i. (8) 

i 

We then proceed to compute the low lying states of 
H e f f using the conventional DMRG again. For this sim- 
ulation we keep m-2 states and use 3 blocks instead 4 to 
form the superblock. 

As expected from the study of S — \ systems, the TS- 
DMRG is more accurate in the highly frustrated regime 
than in the unfrustrated case, pi is relatively large in 
the unfrustrated regime because of the relatively large 
value of the interchain coupling. The same simulation 
with J± = 0.2 leads to an improvment of factor 10. It 
is worth noting that the superblock size in this step is 
mij, we are able to reach ra-i = 100 on a workstation, this 
remains modest with respect to what can be achieved on 
today's supercomputers. For the multi-chain DMRG su- 
perblock sizes of about 100 times larger are accessible^, 
this means that it is possible to reach 7712 ~ 500 on a su- 
percomputer. This would increase the current accuracy 
of the TSDMRG by two or more orders of magnitude. 
This shows the great potential of the TSDMRG. These 
possibilities are under exploration. 



Si = 0^^3,3,0, (4) 
the intra-chain Hamiltonian is likelywise updated 

h = 0\h Bl S 3 ^+h B3 5 hV )0. (5) 

In these equations, we have adopted the usual convention 
that the different blocks of the superblock are labelled 
1-2-3-4. For PBC, blocks 2 and 4 are made of a 
single site and blocks 1 and 3 are the largest blocks. In 
Eq. , it is supposed that the spin to update is in block 
1. In Eq. J5J, liBi represents the internal Hamiltonian of 
block i. The first step ends with the updating of these 
operators. 

Each chain may now be viewed as a super 'spin' with 
additional internal degrees of freedom due to the different 
sites. A chain I is described by its 'spin' value S; = 
(Su,i = l,...,L) and its internal Hamiltonian hi. hi is 



III. RESULTS FOR S = 1 

The case of spin i has been extensively studied in 
Refi 22 i 25 . We were able to reach lattice sizes of up to 
64 x 65 and show that as seen in QMC simulations, the 
system is ordered in absence of frustration for small J±. 
But when Jd ^ 0, we have shown that in the vicinity 
of Jd/J± ~ 0.5 the system is made of weakly-coupled 
chains even when Jj_ and Jd are not small. This find- 
ing of the TSDMRG was checked using ED on small 
systems? 7 . More careful simulations at the vicinity of 
the point Jd/J± = 0.5 revealed that the first neighbors 
interchain correlation, i.e. the transverse bond strength 
is equal to zero at the maximally frustrated point. The 
non-zero correlations, starting from the second neighbor 
decay exponentially. These results lead us to conclude 
that the maximally frustrated point is a quantum crit- 
ical point (QCP) between the two magnetic states (the 
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FIG. 2: Ground-state energy per site as function of Jd for 
J± = 0.4, S = 1 and L = 12. 



second magnetic state is stable when Jd > 0.5JjJ. A 
possible argument against this conclusion is that at the 
maximally frustrated point, the system could be insta- 
ble against higher orders terms such as a ring exchange 
termii In the case of a two-leg ladder this term seems to 
lead to a dimerized state. There are however some strong 
indications that the dimerized state does not exist in this 
model as discussed in Refi^S. The mechanism to avoid 
frustration is to divide the system into chains in which 
the frustrated bonds are severed. We will now study the 
extension of this mechanism to S — 1. We perform the 
same analysis as for S = | for L = 6, 8, 10 and 12. For 
the first three values of L, ED is performed to obtain the 
rri2 lowest eigenvalues and the corresponding eigenstates 
of the chain. For L = 12 the DMRG was used, we kept 
mi = 243 states, i.e, one DMRG iteration was done from 
the L = 10 exact result. The truncation error was about 
7 x 10 -7 . The truncation error is relatively large because 
we used periodic boundary conditions. These rri2 low- 
lying states were then used to generate the 2D lattices, 
i.e, 6 x 7, 8 x 9, 10 x 11 and 12 x 13 respectively. 



A. Ground-state energies 

The curve of E G (J d ) for S = 1, shown in Fig.© is 
similar to that of S — i. Starting from Jd = 0, Eq 
increases until it reaches a maximum at J r d aax . It then 
decreases when Jd is further increased. The position of 
the maximum depends slightly on L and seems to con- 
verge to 0.5Jj_ in the thermodynamic limit. Ec(J d nax ) 
is very close to 1.53, the energy of decoupled chains, but 
always remains slightly lower. Thus as for S — ^ the 
chains are very weakly bound, even though the bare in- 
teractions ( Jj_ = 0.4 and J d = 0.22) are not small. J r d aax 
depends on L as shown in Fig.© and as in the case of 
S = it extrapolates to 0.5 J± in the thermodynamic 
limit. 

Eg (I) shown in Fig.©, where I is the number of chains 
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FIG. 3: Maxima of the ground-state energy as function of L 
for J x = 0.4, S = 1. 
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FIG. 4: Ground-state energy per chain as function of the 
chain number for J± = 0.4, S = 1, L = 10 for Jd — (filled 
circles), and Jd = J™ ax (open circles). 



is dramatically different when J d is far or close to J™ aa: . 
Far from J™ aa: , one of the two magnetic phases is highly 
favored. Starting from an isolated chain, magnetic en- 
ergy can be gained by increasing I , leading to the ordered 
state. The situation is different when Jd — J™ a2; ; neither 
of the magnetic states is favored. At this point, magnetic 
energy cannot be gained and Eq(1) is nearly independent 
of I, the system remains disordered as we will see below 
from the analysis of the correlation functions. 

B. First neighbor correlation 

The transverse first neighbor spin-spin correlation 
taken in the middle of the lattice 

C\ = (Sl/2X/2+iSl/2X/2+2), (9) 

shown in Fig.© is also reminiscent of the S = \ case. C\ 
vanishes linearly at J d = J®- J d 1S slightly different from 
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FIG. 5: First neighbor spin-spin correlation as function of Jd 
for J ± = 0.4, S = 1, L = 12. 



FIG. 7: Transverse spin-spin correlation as function of the 
distance for L = 10, S = 1, J_l = 0.4, Jd = (filled circles) 
and J d = JT ax (open circles). 
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extrapolation of C\ for S = 5 leads to a small negative 
value as shown in Fig. H20|) below. In this case, the order 
parameter is too small to be obtained from an extrap- 
olation from relatively small systems; it is necessary to 
go to larger systems as those studied in Ref.( 2 -) in order 
to extrapolate to the correct thermodynamic limit. The 
extrapolated value of Ci does not however lead to the cor- 
rect value of the magnetization since it is obtained from 
a system with a fixed L. A better estimation is given 
by the finite size analysis of the end-to-center spin-spin 
correlation 



FIG. 6: Minima of the first neighbors spin-spin correlation as 
function of L for J± — 0.4, S = 1. 



J d nax - This small difference is due to numerical error; 
this conclusion is supported by the more accurate results 
obtained for small S — \ systems in Ref.(— ) where J° 
and J™ ax are equal. The extrapolated J° (Fig.© as 
L — > 00 is also in the vicinity of 0.5 Jj_. 



/2,L/2+iSl/2,L + 1) 



(11) 



shown in Fig.©. Cl for L — > 00 is roughly 0.06 consis- 
tent again with the existence of the long-range order. 

In the vicinity of J™ aa: , C; decays exponentially as seen 
in Fig. 10. For 1 = 5, his value is already four order of 
magnitude smaller than in the case Jd = 0. This is con- 
sistent the nearly disconnected chain behavior observed 
for Eq at this point. 



C. Long-distance correlations 

The transverse spin-spin correlation function, 

Cl = (Sl/2,L/2+iSl/2,L/2+i) 



D. 



spin gaps 



(10) 



is shown in Fig.Q for a L — 10 system and J± = 0.4. 
When Jd = 0, C; extrapolates to a finite value in the 
thermodynamic limit. This result is important because 
it shows the non-perturbative nature of the TSDMRG. 
Starting from an isolated chain which is disordered, the 
TSDMRG can reach the ordered phase. The ordered 
phase can be easily reached for larger spin than for S = 5 
where quantum fluctuations are more important. The 



The variation of the spin gap A with I (for a fixed L), 
L, and Jd the are also consistent with the above findings. 
A(J<j) for L = 10 system is shown in Fig.©. A(J d = 0) 
is about 0.05 (in this regime the gap is zero in the thermo- 
dynamic limit as we will see below); A remains relatively 
flat as Jd is increased until it reaches the vicinity of J d nax . 
Near J™ a2: , A(Jd) first sharply increases and reaches the 
finite size gap of an isolated chain. As Jd > J d nax , A( Jd) 
first sharply decreases and then becomes nearly constant 
at about 0.05. A (I) (Fig.^DJ) is reminiscent of Eq(1); in 
the unfrustrated case, the chains are effectively coupled. 
A(l) rapidly decreases from about 0.53 to 0.05 as I is var- 
ied from 1 to 11. At Jd = J d nax however, A(Z) is nearly 
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FIG. 8: Center-to-end spin-spin correlation as function of L 
for S = 1, J± = 0.4 and J d = 0. 



FIG. 10: Gap as function of the number of chains for L — 10, 
S = 1, Jl = 0.4, J d = (filled circles) and J d = JJ 1 ** (open 
circles) . 




FIG. 9: Gap as function of Jd for L — 10, S — 1; the dotted 
line represent the single chain gap. 
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FIG. 11: Gap as function of L for S = 1, Ji = and Jd = 
(single chain, filled circles), J± — 0.4, Jd — (filled squares) 
and J d = J™ ax (open circles) . 



independent of I. 

The analysis Cl shows that for Jd — 0, the system 
is ordered. We thus expect that A(L) — > as L — > 
oo. This is seen in Fig.Qll)) where A(J) is shown for 
L = 6,8,10 and 12 systems. The decay faster than ^ 
the extrapolation to leads to a negative value. At Jd — 
jmax on orc j er hand, A(L) remains close to that of an 
isolated chain in all case. The two functions are finite in 
the thermodynamic limit. This result show the dramatic 
difference between S = \ and 5 = 1 systems. For S = A, 
an equivalent plot lead to a zero gap2£ at the maximally 
frustrated point. The extrapolated value for J d = 0, A = 
0.4015 agrees well with the current best estimate of the 
Haldane gap A# = 0.4107. The difference is due to the 
relatively short chains, up to L = 16, that were used for 
the extrapolation not to the DMRG that yielded highly 
accurate results for each size studied. Noting that the 
spin-spin correlations in the transverse direction have a 
very short range, the difference between the extrapolated 
values of J d — and J d — J d nax appears to be relatively 
large. We believe that this difference could be inferred 
from the fact that the extrapolation from 2D systems 



were done with lattice sizes up to L — 12 only. 

E. Spin-spin correlations on large systems 

So far, in the study of S = 1 systems, we have fixed 
Jl = 0.4. This choice was motivated by the presence of 
Ah = 0.4107 in an isolated chain. We initialy felt that 
it was necessary to choose a large enough J± so that the 
chains will effectively be coupled when the perturbation 
is turned on and this will lead to sizable correlation in 
the thermodynamic limit. But, this choice limited us to 
relatively small lattices, L < 12. This is because when 
Jl is large, the condition 5E/J± ^ 1 is hard to fulfill 
for larger L. For instance for L = 16, SE/J± = 5.65 for 
m 2 = 64, this prevented us to study L = 16 lattices. But 
in the course of this work, we find that even smaller values 
of Ji can lead to detectable values of in the unfrustrated 
regime C; as I — > oo. For smaller Ji, we can actually 
reach larger L. We wish to present in this part our results 
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FIG. 12: Longitudinal spin-spin correlation as function of the 
distance for L = 24, S = 1, J± = 0, Jd = (filled circles), 
J± = 0.2, J d = (filled squares), J d = 0.2, J d = 0.102 (open 
squares) . 



FIG. 14: Magnetic structure factor for L = 24, S — 1, Jj_ 
0.2, J d = 0. 
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FIG. 13: Transverse spin-spin correlation as function of the 
distance for L — 24, S = 1, Jj_ = 0, Jd = (filled circles), 
J± = 0.2, J d = (filled squares), J d = 0.2, J d = 0.102 (open 
squares) . 



for Ji = 0.2 and L = 24. These results will add strength 
to those of J± = 0.4 presented above. 

In Fig. l|12l) . we show the longitudinal correlation, 



/2+1,L/2+iSl/2+1,L/2+1/i 



(12) 



taken in the middle chain. For J± = 0.2 and Jd = 0, 
Ci clearly extrapolate to a finite value as expected. But 
for J± = 0.2 and Jd = 0.102, C\ is nearly identical to 
spin-spin correlation on an isolated chain. For the trans- 
verse correlation C; shown in Fig. (|13fl . we see again the 
dramatic difference between the unfrustrated and highly 
frustrated cases. In the first case, Ci goes to a finite 
value when I — > oo. But for the highly frustrated case Ci 
decays exponentially. 

Another picture of this dramatic difference is given by 
the magnetic structure factor S{Q = (q x ,q y )) shown in 
Fig. (|14ll5|l . In the magnetic phase, S(Q) is dominated 
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FIG. 15: Magnetic structure factor for L 
0.2, Jd = 0.102. 



24, S = 1, Ji 



by a sharp peak at Q = (tt, it) indicative of Neel order. 
In the disordered phase, S(Q) is nearly flat except for a 
small range along Q = q y ) which retains the signature 
of short-range in chain correlations. 



F. Conclusion 

In this section, we have presented comprehensive re- 
sults on S = 1 coupled chains. These results show some 
analogy with those of S — | published in Ref.( 22 ). Start- 
ing from the unfrustrated system for which Jd = 0, the 
ground state for S = 1 is ordered, as expected from the 
DLS theorem. While for S = h it was necessary to sim- 
ulate lattices of up to L = 64 25 in order to see the ex- 
trapolation of Ci to a finite value, for 5 = 1, relatively 
short sizes (L = 12) were enough. This is indeed because 
quantum fluctuations are less important in a S = 1 sys- 
tem, i.e., the order parameter is larger. This enables it to 
be computed more easier. By comparison, the same ex- 
trapolation done for S — \ will lead to a negative value. 



When the frustration Jd is turned on and reaches the 
value J™ ax , a point where Eq is maximum, Cj decays ex- 
ponentially, A takes a value very close to that of a pure 
ID system. These results imply that at J™ ax ) the trans- 
verse interactions are irrelevant. For 5 = ^, we identified 
this state as spin version of an SLL. In that case, this slid- 
ing phase will probably be confined at the critical point 
where the two competing magnetic states Q = (tt, tt) and 
Q = (tt, 0) neutralizes each other. However, we cannot 
completely rule out a small finite extension of the sliding 
phase or even totally exclude the emergence of a relevant 
interaction at lower energies which eventually drives the 
system to a dimcrizcd phase 13 . It is obvious that, though 
both 5 = | and 5=1 systems are made of nearly discon- 
nected chains at J™ aa: , the conclusions must be different 
because of the presence of the Haldane gap Ah in the 
5 = 1 chain. The existence of Ah restricts the possible 
phases that may arise in the vicinity of J d nax . The first 
crucial difference is that the disordered state that exist 
in an 5 — 1 system has a gap in its excitation spectrum 
as seen in Fig. (jl 1|> ; any eventual residual interaction will 
be wiped out by this gap which means the emergence of 
new phases at low energies is not favorable for an 5 = 1 
system. This disordered phase has probably a finite ex- 
tension which is roughly oc Ah- 



IV. RESULTS FOR S = \ TO S = 4 

In the study of various 5, we will not do the same 
extensive calculation seen in the preceeding section with 
5 = 1. We will simply fix L and analyze the behavior 
of the system as function of Jd and I. As we will see 
below, quantum fluctuations are small for 5 > 1, the 
study of relatively small systems is enough to get the 
correct picture in the thermodynamic limit. We studied 
a lattice with L x (L + 1) = 10 x 11 for 5 = \ to 5 = 4. 
J± was set to 0.4 so that it is larger than the finite size 
gap in spin half-integer systems and larger or close to the 
Haldane gaps in spin integer systems. Jd is varied from 
to 0.4. 



A. Ground-state energies 



The ground-state energy is shown in Fig. i|ltj|) for 5 

7 



2 

to ^. Simulations were also done for 5 = 4 but they 
did not converge for certain values of Jd- We infer this 
failure to the large degeneracy of the renormalized single 
chain Hamiltonian hi for large 5. At Jd = 0, the curves 
approach the classical value Eq/S 2 = — 1.4 quite rapidly. 
For 5 = |, we find E G /S(S = 1) = -1.093. But if we use 
the 1/5 2 normalization, we get Eq/S 2 = —1.405. Hence 
if both finite size effects and the correct normalization 
are taken into account, 5 = | is already in the classical 
limit. 

Larger spin systems are found to present the same fea- 
tures displayed by spin 1/2 systems. The ground-state 




FIG. 16: Ground state energy as function of Jd and 5 — | 
(filled circles), S = 1 (open circles), S — | (filled squares), 
5 = 2 (open squares), S — | (filled diamonds), 5 = 3 (open 
diamonds), 5 = |, (filled triangles). 




FIG. 17: Ground state energy at Jd = as function of I and 
5=| (circles), 5=1 (squares), S = | (diamonds), 5 = 2 
(triangles) 



-5 



-6<K 



? -7, 



-8"- 



-10 



10 



12 



FIG. 18: Ground state energy at Jd = as function of I and 
5 = i (circles), 5=1 (squares), 5 = § (diamonds), 5 = 2 
(triangles) 
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FIG. 19: First neighbor correlation as function of Jd and 
S = | (circles), S = 1 (squares), S = | (diamonds), S — 2 
(triangles) S 1 = § (filled diamonds), 5 1 = 3 (open diamonds), 
S = I, (filled triangles). 



energy, Eq, shown in Fig. (|16|> increases as Jd increases 
until the maximally frustrated point where Eq of the 
two-dimensional system becomes very close to that of 
disconnected chains. From this point it decreases when 
Jd is further increased. This may be interpreted as fol- 
lows: starting from the Neel state with Q = (tt,tt) for 
Jd = 0, the system tends to lose energy under the ac- 
tion of Jd which progressively destroys the Neel order 
until the maximally frustrated point. Beyond this point, 
Jd becomes dominant and the systems enters the Neel 
Q = (n, 0) phase. The position of this maximum de- 
creases slowly with increasing S. This indicates that in 
addition to the effect of OBC that shifts J™ ax towards 
higher values, there are intrinsic finite size effects. All 
systems evolve regularly towards the S — * oo limit. 

The curve of Eq appears to change structure as S 
increases. At low S, a well-rounded maximum is ob- 
served. We were able to fit all the points of the curve to 
a quadratic function. But for large S, this became im- 
possible. The maximum has nearly become a cusp as for 
S — ► oo. This cusp is at the intersection of two straight 
lines E Gl = -1 - J± + 2J d and E G2 = -1 + J± - 2J d 
which are the ground-state energies, respectively, below 
and above the transition point J d = 0.5Jl- 

Ec(l) shown in Fig. (|17ll8fl displays the features seen 
for S = 1 . It decreases when I increases in the weak frus- 
tration regime. It remains nearly constant in the vicinity 
of the maximally frustrated point. 



B. First neighbor correlation 

The tendency to the severing of the chains is more 
clearly seen in the transverse bond strength 



Ci = (S z 5 , 6 Sl 7 ), (13) 



0.25 




FIG. 20: Spin-spin correlation at Jd = as function of I 
and S = | (filled circles), S = 1 (open circles), S = | (filled 
squares), S = 2 (open squares), 5=1 (filled diamonds), 
S = 3 (open diamonds), S = |, (filled triangles), S — 4 (open 
triangles). 



shown in Fig. l|19|) . In all cases, C\ decreases from its 
value at J d — and seems to vanish at J° (we were 
able in all cases to reach values of C\ which are equal or 
less than the numerical accuracy of in our simulations). 
From this point it increases. There is a small difference 
between the position of J° for different values of S as 
found for Eg- The curves of C\ suggest that for all S, the 
mechanism to avoid frustration is identical: the systems 
relaxe to nearly disconnected chains. These curves also 
show the influence of quantum fluctuations for small S. 
This is seen in the decay of C\ as soon as Jd ^ 0. For 
larger S, C\ remains nearly constant until J d J°. 

C. Long-distance correlations 

Since our starting point for 2D systems is disconnected 
chains, it is important to show, as for the spin 1/2 case 
studied previously, that the TSDMRG is able to reach 
the ordered phase. One possible way to look at the ap- 
pearance of the ordered state is to look at the decay of 
the transverse correlation function, 

Cl - (Si e Sl 5+l ), (14) 

in the Neel phase for Jd = 0. We found that for all values 
of S except S — \, as shown in Fig. (|20|) . the transverse 
correlation function extrapolate to finite values. C; ex- 
trapolates to a negative value for S = \. In that case, 
quantum fluctuations are so strong that it is necessary to 
go to larger values of L as done in Ref.(— ). 

At the maximally frustrated point, we also observe 
similar exponential decay of Ci as for S = \. This de- 
cay is less faster with increasing S. Indeed in the limit 
S — > oo, the transition is of first order. The chains are 
disconnected in this classical limit and C/ is exactly equal 
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FIG. 21: Transverse spin-spin correlation at Jd = 0.22 as 
function of / and S = | (circles), S = 1 (squares), S = § 
(diamonds), S 1 = 2 (triangles), S — § (filled diamonds), S 1 = 
3 (open diamonds), S = |, (filled triangles), S 1 = 4 (open 
triangles). 



to zero. However for any small deviation from the tran- 
sition point, the system falls in one of the ordered states. 
This point is virtually impossible to find exactly numeri- 
cally. However, for smaller values of S the critical region 
is larger and even if we miss the exact transition point, 
this behavior will nevertheless be observed as far as we 
close enough to the QCP. 



D. Spin gaps 

The curves of A(Jd) in Fig. (|22|l for different values of 
S have typically a peak in the at Jd — J r d aax . This peak 
is very narrow, except for S = \ where quantum fluctu- 
ation effects lead to a broader peak. As expected from 
the behavior of C\ , this peak is sharper with increasing 
S. A(J™ aa: ) is nearly equal to the finite size gap of an 
isolated chain which is represented by a flat line in each 
case. We were unable to reach the ID gap for S > f . 
This is probably due to the narrowness of the critical 
region which makes it difficult to see the nearly discon- 
nected chain regime. One can easily fall in one of the 
ordered regime leading to a relative slower decay of Ci 
which manifests itself to a smaller finite size gap. 



E. Conclusion 

In this section, we presented results for L = 10 and 
J± = 0.4 with S varying from ^ to 4. In agreement 
with the DLS theorem, we found long-range order for 
all S greater or equal to 1 in the unfrustrated case. As 
Jd is turned on, the long-range order is destroyed. An 
interesting question is the nature of this disordered state 
as function of S. Before addressing this question, we will 
first review how frustration works in ID 38 . 




FIG. 22: Spin gap as function of Jd and S — \ (circles), S = 1 
(squares), S = | (diamonds), S — 2 (triangles) 



For the frustrated S = | chain, there is a transition 
to a dimerized phase at J^ c = 0.242J 1; where J 2c is the 
value of next-nearest neighbor coupling at the critical 
point. At this point a gap opens exponentially and grows 
with J 2 . At J2 = 0.5Ji, the system is perfectly dimerized 
and shows incommensurate correlation above this point 
(disordered point). At about J2 = 0.52Ji a two-peak 
structure appears in the structure factor (Lifshitz point). 
DMRG simulations for the S = | chain shows a similar 
behavior. This suggests that it is generic to half-integer 
spin systems. 

Integer spin chains are already gapped in the absence 
of the frustration term J2. For S = 1, the transition 
to a dimerized phase is absent, but numerical simula- 
tions show the presence of the disordered and the Lif- 
shitz points. In addition there is a first order transition 
at J2 = 0.75Ji from a phase with a single string order 
to a phase with a double string order. At this point the 
chain splits into two chains. These special points are 
also observed in the S = 2 chain except that the order 
parameter for the first order transition is still unknown. 

The results presented in the preceding sections show 
that the mechanism to case frustration works differently 
in 2D systems. This mechanism is the same for all val- 
ues of S. The system spontaneously severs the frustrated 
bond at the maximally frustrated point. The similarity 
for all 5* of this mechanism stems from the fact that if 
the transverse coupling is large enough, all the 2D sys- 
tems are ordered for half-integer as well as for odd integer 
systems. Frustration is a competition between two mag- 
netic ground phases and we have shown that for coupled 
chain systems, the best way to avoid frustration is to re- 
lax into nearly independent chains. It is clear that such 
a mechanism will be independent on the value of the spin 
as found in our numerical study. The consequences are 
nevertheless different for half-odd integer and for integer 
spin systems. 

For all spin half-odd integer systems, like for the spin 
1/2 studied more extensively in Ref.(— ), there is a second 
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order phase transition between the two magnetic states 
at J d = J™ ax . At the critical point, the system is dis- 
ordered. The transverse correlation decay exponentially 
while at long distances, the longitudinal one behave like 
those of independent chains. Hence at the critical point, 
the spin rotational symmetry of the Hamiltonian is re- 
stored. As for spin 1/2, there might be a residual in- 
teraction which can drive the system eventually to a SP 
phase. But, previous numerical studies on 2D systems^ 
and on three-leg ladders point to an absence of a dimcr- 
ized phase in this region for S = ^. This is expected 
to be valid for all half-odd S. We would like to stress 
that dimerization is not the driving mechanism mech- 
anism in the formation of the disordered state. We are 
indeed aware of earlier ED results^ in which an enhance- 
ment of the SP susceptibility was observed in the regime 
J 2 w Ji/2. We believe in the light of our results that this 
merely the consequence of the severing of the chains in 
one of the two directions of the square lattice. The SP 
signal is expected to be larger in ID where it has a power 
law decay than in 2D when the spins are locked into Neel 
order in the unfrustrated regime. 

For integer spins, there is an intermediate phase be- 
tween the two magnetic states. When \ J± — 2Jd\ <C A, 
where A is the single chain spin gap, the transverse cou- 
pling are irrelevant. The maximally frustrated point is 
the equivalent of the disordered point seen in ID. In this 
regime of couplings, the system is an assembly of nearly 
decoupled chains. In the case of integer spins, even if 
there is a residual interaction at the maximally frustrated 
point, this interaction is necessary irrelevant because of 
the presence of A. Integer spin systems are thus radically 
different from half-odd integer systems. 



V. CONCLUSION 

In this paper, we used the TSDMRG to study coupled 
spin chains with S varying from i to 4. This study illus- 
trates the power of the TSDMRG method, where using 
a modest computer effort we were able to study the un- 
frustrated regime and find long-range magnetic order, in 
agreement with the DLS theorem and Monte Carlo stud- 
ies. We obtained good accuracy in the highly frustrated 
regime of the model. The study of this region has so far 
resisted to other numerical methods. 

We showed that in order to avoid frustration, all spin 
systems tend to sever the frustrated bonds. The severing 
of the transverse bonds is a large effect which is seen 
in various physical quantities. The strong frustration 
regime is dominated by ID physics, topological effects 
become important as predicted in Ref.(i&2i). However, 
we did not find any qualitative difference between odd 
and even integer spin systems as predicted in Ref.(2*2I). 
It could be due to the fact that in the highly frustrated 
regime the 2D systems tend to relax into nearly indepen- 
dent ID sytems where tolological effects are identical for 
odd and even integer spins. It could also be related to 
the anisotropy of the model studied. 
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